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We study the nonequilibrium dynamics of an electronic model with competing spin-density-wave 
and unconventional superconductivity in the context of iron pnictides. Focusing on the collisionless 
regime, we find that magnetic and superconducting order parameters may coexist dynamically 
after a sudden quench, even though the equilibrium thermodynamic state supports only one order 
parameter. We consider various initial conditions concomitant with the phase diagram and in a 
certain regime identify different oscillatory amplitude modes with incommensurate frequencies for 
magnetic and superconducting responses. At the technical level we solve the equations of motion for 
the electronic Green’s functions and self-consistency conditions by reducing the problem to a closed 
set of Bloch equations in a pseudospin representation. For certain quench scenarios the nonadiabatic 
dynamics of the pairing amplitude is completely integrable and in principle can be found exactly. 
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I. INTRODUCTION 

Conventional superfluids and superconductors host 
various collective oscillations. The best studied examples 
include the phase mode of the order parameter (OP), the 
so-called Anderson-Bogoliubov mode,li^^ and the ampli¬ 
tude oscillations in the magnitude of the superco ndu ct- 
ing gap, the so-called longitudinal Schmid mode.^^ In 
charged superfluids the coupled oscillations in the phase 
of the order parameter and the electric field appear be¬ 
cause of gauge invariance. Physically, this mode corre¬ 
sponds to the balanced oscillations between the normal 
current and supercurrent, and in the lite ratu re it is called 
the transverse Carlson-Goldman mode.l^ Early works 
on superconducting modes were comprehensively sum¬ 
marized by Artemenko and Volkov,^ and Kulik, Entin- 
Wohlmant, and Orbach,!^ including studies of disorder 
scattering effects and charge imbalance on the dispersion 
and attenuation of collective oscillations. 

In multicomponent systems or superconductors with 
unconventional symmetry of the OP, the plethora of col¬ 
lective effects is even richer.l^®^] multiband supercon¬ 
ductors such as MgB2 the oscillations of the phase dif¬ 
ference of OPs between the two bands is charge neu¬ 
tral, in contrast to phase average plasma oscillations. A 
phase difference low frequency Leggett mod^^^ is an in¬ 
gap weakly damped excitation observed in the Raman 
response of MgB2.^ It is natural to look for Leggett- 
like and Carlson-Goldman modes i n mu ltiband and iron- 
pnictide superconductors (FeSCs).^^^^ Normally, these 
modes are overdamped with frequencies well within the 
quasiparticle continuum. There are, however, important 
exceptions to this rule. A typical setting for this scenario 
is the change in the OP symmetry controlled by external 
parameters. In many cases such a transformation pro¬ 


ceeds via an intermediate phase with broken time reversal 
symmetry (TRS).fi 2 HIIl Soft Leggett-like modes are found 
at the boundaries of the intermediate lower symmetry 
phase. A transformation of this kind was very recently 
induced by pressure in KFe2As2,^ and TRS breaking 
along with Leggett-like modes await experimental detec¬ 
tion. 

Different kinds of c ollecti ve excitations are the 
Bardasis-Schrieffer modes.l^SlSni These modes are carried 
by Cooper pairs accelerated to higher angular momen¬ 
tum states and manifest as in-gap excitons. Since all 
but the s-wave channels are charge neutral, Bardasis- 
Schrieffer modes remain low energy in-gap excitations 
even in the presence of Coulomb repulsion. As photons 
transfer the ang ular m omentum to Cooper pairs, Ra¬ 
man spectroscop \ l^d 32 | jg ideally suited to probe Bardasis- 
Schrieffer modesP^lf^ 

Interestingly, Raman spectroscopy was originally sug¬ 
gested as a tool to detect amplitude Higgs modes, 
whereas they were indirectly observed owing to cou¬ 
pling to intermediary collective excitations.^SHini More 
recent Raman,!^ terahertz pump-probe spectroscopy,^^ 
and combined tunneling and optical conductivity mea¬ 
surement^^ provide unambiguous direct tests of mas¬ 
sive Higgs modes in superconductors. In a parallel vein, 
the coherent amplitude mode has been observed in the 
strongly interacting superfluid phases of bosonic ultra¬ 
cold atoms in optical lattices b y mea ns of Bragg spec¬ 
troscopy and lattice modulation. b^ l ^^ l All these findings 
have stimulated many theoretical efforts, (see a recent 
review article in Ref. | 46 ) and references therein). 

Since the pioneering work by Volkov and Kogan,!^ per¬ 
sistent oscillations of the superconducting OP have been 
predicted to appear in a response to a fast nonadiabatic 
perturbation (quench) in the collisionless regime.l^^t^ 
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The mode frequency is determined by the supercon¬ 
ducting gap whereas oscillations are superimposed with 
a slow power-law decay. In contrast, theoretical stud¬ 
ies of nonequilibrium dynamics after ultrafast excitation 
in complex superconducting systems hosting coexisting 
OPs are in their early stages, with only a few recent re- 
s^ltsIUHIIl The main thrust of this paper is to provide 
a detailed description of the coupled dynamics of am¬ 
plitude modes in the context of FeSC systems. Broadly 
formulated, our theory may shed light on the hotly de¬ 
bated issue of the structure of the OP and the closely 
related question of the competition between magnetism 
and superconductivity in FeSCs as seen out of equilib¬ 
rium. Our motivation comes from recent ultrafast mea¬ 
surements that reveal a dynamic coexistence and inter¬ 
play of multiple order parameters in various strongly cor¬ 
related materials 

This paper is organized as follows. In Sec.|n]we adopt 
the simplest model of iron pnictides, where the dynamics 
of competing orders is already found to display a nontriv¬ 
ial character. We briefly discuss the ground state prop¬ 
erties of this model and derive the equations of motion 
for the Green’s functions in the pseudospin representa¬ 
tion. In Sec. |III| we numerically integrate these coupled 
equations together with self-consistency constraints, and 
discuss the emergent dynamic coexistence of supercon¬ 
ductivity and magnetism. Section]^ is devoted to the 
analysis of a special case when the dynamics of the order 
parameters is integrable. We summarize our results in 
Sec. |V] and place our work into the perspective of future 
studies. 


II. MODEL 

To study th- dynamical interplay between supercon¬ 
ductivity and spin-density wave order, we use the mini¬ 
mal two-band model previously introduced in the context 
of iron-pnictide superconductors. Following Refs, 
we consider the Hamiltonian 

H = H^ + H^ + H^. (I) 

The first term accounts for the electronlike and holelike 
two-dimensional (2D) fermionic bands, 

Hq = |CkcC]{jQ,Cka + Ck'//k'Q;/k'a| i (2) 

k 

where /k/a, /k'a (k^ = k—Q) are the creation and annihi¬ 
lation operators for the fermions with a spin projection a 
near an electron pocket Q = (0, tt) of the two-dimensional 
Brillouin zone with dispersion 

^2 

fk/ = y - M/, (3) 

the chemical potential /i^, and we set the electron’s mass 
to one. Similarly, Cka describe the fermions near 


the F = (0, 0) point with the hole band with a chemical 
potential [ic and dispersion 

^kc = Me (4) 

The second term in Eq. Q accounts for the supercon¬ 
ducting pairing. Within the mean-field approximation 
we have 

=2 {^a/3Ckc<clj^^ -I- + h.C.j , ( 5 ) 

k 

where are the superconducting order parameters de¬ 
fined for each band, 

= gsc'^{icry)ai3{iay)\s{f-k^fi,s), 

( 6 ) 

= 9sc'^{icry)ai3{iay)\s{c-UyCks), 

k 

and gsc > 0 is the superconducting coupling. Finally, the 
last term in Eq. Q describes the onset of the commen¬ 
surate spin-density-wave (SDW) order, which within the 
mean-field approximation is described by 

Hm = 2 ■[/ka'^k/S + Ckct/k/?! + h.C., (7) 

k 

where the SDW order parameter is determined self- 
consistently via 

( 8 ) 

p 

and Mm > 0 is the corresponding coupling constant. 
In what follows, without loss of generality, we assume 
that (i) there is no mismatch between the hole and 
electron Fermi surface, Mc = M/i (ii) Wa/3 = 
and (iii) the superconductivity is of type, so that 
= Ac,/(icr^)a/3 and 

Ac = -A/ = A. (9) 

Below, we briefly review the ground state properties of 
the mean-field model Q for the superconducting pair¬ 
ing. 

A. Correlation functions and pseudospins 

In order to obtain the equations of motion, which de¬ 
termine the nonadiabatic dynamics for superconducting 
and magnetic order parameters, we introduce the four- 
component spinor 

^ka ~ ( ''ka’ C-ka, /kc,, /-ka ^ • (10) 

The corresponding real-time correlation functions are 
Ga;3(k;tl,t2) =-*(F{^ka(tl)4'^(^2)}). (11) 




3 


We further consider normal-G and anomalous--F' propa¬ 
gators for c-fermions, 


G°^(k;ti,t2) = -i(T’{cka(G)4^(t2)}), 

= -*(r{cka(G)C-k/3(t2)}), 
Flpiy-MM) = -*(T’{cLk„(ii)c^(i2)}), 
G°^(k;ti,t2) = -i(r{cLkc(G)c-k/3(i2)}), 

and analogously for /-fermions. In addition, we also con¬ 
sider the mixed correlators 

G-^(k,q;G,t2) = -*(T{ck„(G)/^^(<2)}), 
G-^(k;G,t2) = -^{f{fUtl)<ipit2)}). 


The correlation functions above depend on ti and t 2 . 
However, the magnetic and superconducting order pa¬ 
rameters are determined at ti = ^2 and therefore will 
depend on t = (G -|- <2)/2 only. Accordingly, we intro¬ 
duce the following pseudospin variables P 


defined according to 

L,(k,t) = -1 ^ {G=^(k;t) + GL(k;t)}, 

= \ sign(a)|G^5(k;t)-f G^g(k;t)|, 

L^(k,t) = -^ ^ sign(a) {G^„(k;t)-f G;^„(k;t)} . 
a=td 

(16) 

Equations of motion for the pseudospins can be ob¬ 
tained from the equations of motion for the fermionic 
operators. Using the Heisenberg representation 

Cka(i) = (17) 

we have 
d 

f^Cka = ^c(k)ckQ -|- ka ^ ^ nZct/j/k/Ji 

d -c "" 

*^'^k;/3 “ ~?c(k)c^ — ^ ,/kQnT-a/3, 


if/(k,t) = /f/(k,t) - iKy{\^,t) = *E/^(k,t), 

7^7 (k,f) = - *i^}(k,t) = *E/^(k,t), 

7s:+(k,t) = fF4(k,t), A:+(k,t) = fF^^(k,t), (14) 

Y. Gaa(k:i)- 

a=T.t 


where we used /c(k) = /c(—k). Similarly, the equations 
for the /-operators are 


*^/ka = C/(k)/ka + ^aaf-Wa + "mapCkp, 


dt 

d 


ai/k/3 — 0(k)/k,a ^p/3f-k/3 


(19) 


dt 


Cka’^a/: 


Similarly, we introduce the additional pseudospins <S'(k, t) 
and iV(k, t), which are defined by the mixed averages 
from Eq. (p^: 


Here A denotes the complex conjugate of A and a = —a. 
Erom these equations we derive the equations of motion 
for the correlators above. 


S,{k,t)='- Y {GSj(k;t)+GS,(k;t)}, 
a=td 

Sy(k,t) = -^ Y sign(a) |G((7(k;t)-hG(75(k;t)|, 
a=td 

5,(k,t) = ^ ^ sign(a) {G™Jk;t) + G™ (k;t)} , 
a=td 

A,(k,t) = i 5] {G™5(k;t)-G-s(k;t)}, 
a=td 

NyiKt) = ^ Y sign(a) {G((7(k;t) - G((7(k;t)| , 
a=td 

Nz(k,t) = ^ Y sign(a) |G“„(k;f) - G((7(k;t)| . 
a=td 

(15) 


B. Ground state 

In this section we discuss the ground state proper¬ 
ties of the mean-field Hamiltonian Q. Generally, the 
ground state properties can be derived by analyzing the 
free energy within the Luttinger-Ward’s generating func¬ 
tional method.l^ For convenience we adopt the pseu¬ 
dospin variables. The distribution of the pseudospins 
variables ATkc./ with respect to momentum follows from 
the mean-field theory. For example, for the ^ we find 

^ — /k/,c (Ckc,/ 4" 4" |A’k+E’k-|^ 

2p;2 JgkTl + ■ 

( 20 ) 

Here we have introduced the renormalized quasiparticle 
spectrum 


Finally, we will also need pseudospins L(k, t), which are 


4?k± = yei + <5^ + m2+A2±2|dkiym2^ (21) 
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with = (^k/ - ?kc)/2 and function dk = (^k/ + ?kc)/2 
which accounts for the Fermi surface mismatch. The 
anomalous x-components of the K are given by 

2(Bj^|Bk-| + B2JBk+l) ’ 


and the remaining components are zero. Setting in the 
equations above <5k = 0, the pairing amplitude is deter¬ 
mined by the corresponding self-consistency condition 

^ = ~9sc Kpc — 9sc 

p p 



where 

Similarly, pseudospin variables S'k are determined by 
the mixed correlators. In the ground state = 0 and 

qz _ (l^'k+Ii'k- I + — ^k/Ckc + 

El^\E^^\+El_\E^^\ ■ 

The SDW order parameter is determined by 

mz = gm^S^. (25) 

k 

Numerical analysis of the self-consistency equations 
above shows that for the zero Fermi surface mismatch 
the ground state is determined by the ratio of the cor¬ 
responding coupling constants for the magnetic and su¬ 
perconducting orders. When the critical temperature of 
the SDW transition is higher than the superconducting 
critical temperature, Tm > Tc, the system orders mag¬ 
netically, rriz 0, A = 0. If one allows for a finite Fermi 
surface mismatch, then there is a critical value for Jk 
when the superconducting order becomes energetically 
favorable. Furthermore, the thermodynamic phase dia¬ 
gram contains an interm ediate region where two order 
parameters coexist 


C. Equations of motion 

Pseudospin variables happen to be very convenient to 
describe the nonequilibrium dynamics of the magnetic 
and superconducting order parameters. The same tech¬ 
nique has been recently employed to describe the Higgs 
mode in conventional superconductorsWe find that 
the dynamics can be fully accounted for by the five pseu¬ 
dospins, which have three components each. Pseudospins 
Kcj(k, t) describe the dynamics of the electronic degrees 
of freedom of hole and electron bands respectively. The 
corresponding equations of motion are: 

dtKpc = Bpc{t) X Kpc{t) + (m{t) ■ Np{t)) , 

> . X (26) 

dtKpf = Bpf(t) X Kpf{t) — Cz {m(t) ■ Np{t)j , 


FIG. 1: (color online) Superconducting and magnetic order 
parameter dynamics for the metastable initial conditions for 
the choice of coupling constants corresponding to the ratio of 
critical temperatures Tc = 1.75Tm- Both A(t) and m{t) are 
given in units of the Fermi energy £ f - We consider N = 1004 
pseudospins. The bandwidth £a = IOef and the level spacing 
is 5 = £a/N. 


where we introduced the effective magnetic fields 

Bpa = 2{-Al{t),-Ayit),^pa), (27) 


and Ak(t) accounts for the influence of the magnetic or¬ 
dering on the superconducting dynamics. The remaining 
equations of motion are 


dtSp + 2^pNp{t) + 2m(t) x Lp{t) = 0, 

dtLp -I- 2m(t) X Sp{t) = 0, (28) 

dtNp - 2^pSp{t) + 2rh{t) = 0. 


Equations (26) - (28) represent the main result of this 
section. We will analyze these equations numerically in 
Sec. III. In Sec. IV we consider the special case when 
the first two equations (26) decouple from the rest: this 
situation corresponds to quenches of the magnetic cou¬ 
pling constant to zero. In this scenario, the equations 
of motion can be integrated exactly and we prove the 
integrability of this particular case. 


III. DYNAMICAL COEXISTENCE OF 
SUPERCONDUCTIVITY AND SDW ORDER 


In this section we solve the equations of motion (26) 


and (28) numerically. We consider the initial conditions 


corresponding to the metastable state of coexisting mag¬ 
netism and superconductivity, ruz 0 and A y^ 0. 
For the K and S we choose the initial configuration 
corresponding to the metastable state where both rUz 
and A are nonzero, Eqs. (20), (22), and (24). In ad¬ 


dition, for the initial conditions m = rUzCz, we find 
Lk(t = 0) = iVk(t = 0) = 0. 
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We present the results of the numerical integration 
of the equations of motion in Figs. HH In Fig. we 
choose the parameters corresponding to the supercon¬ 
ducting ground state, = l.TSTm. We see that in this 
case magnetization vanishes dynamically, while the pair¬ 
ing amplitude remains finite. 

The results in Fig. [^were obtained for = 0.95Tm- In 
this case we hnd that both magnetization and the pair¬ 
ing amplitude coexist dynamically. We observe that this 
nonequilibrium effect persists for the range of parameters 
corresponding to Tc ~ Tm- 

Finally, in Fig. we show the time evolution of the 
pairing amplitude and magnetization when the initial val¬ 
ues of magnetization and the superconducting energy gap 
are such that Tc = Q.bTm. In this case, we see that the 
pairing amplitude vanishes dynamically, while magneti¬ 
zation remains finite. 

The first and third scenarios are similar to previ¬ 
ously studied cases of collisionless relaxation in a single¬ 
component system. The surviving order parameter, 
which corresponds to a thermodynamically favorable 
state, exhibits oscillatory behavior superimposed with a 
rather slow power-law decay at long times. The physical 
mechanism of relaxation is analogous to collisionless Lan¬ 
dau damping in plasmas. The difference, however, is that 
in a gapped system such relaxation is typically nonexpo¬ 
nential because of the branching singularity in the density 
of states. For example, the superconducting response was 
shown to fall asymptotically as oc cos(2At)/-\/t. 

The second scenario in Fig.j^is special and representa¬ 
tive of the case when both order parameters are of com¬ 
parable strength. Then the initial thermodynamically 
metastable state survives out of equilibrium for extended 
times until the system enters into the collision-dominated 
regime of relaxation. Depending on the choice of param¬ 
eters in the model, both order parameters may oscillate 
with incommensurate frequencies. 

The dynamical effect of coexistence has been recently 
pointed out in the case of multiband superconductors, 
which can be applicable to either MgB 2 or iron-based 
superconductors in the part of the phase diagram with¬ 
out magnetism.l^ In this case mutual dynamics is pri¬ 
marily triggered by the Josephson coupling of pair am¬ 
plitudes between the bands. This is quite different as 
compared to the model which hosts order parameters 
whose physical nature is not the same. Known exam¬ 
ples of the latter kind include the dynamical coexis¬ 
tence of bond-density-wave and d-wave superconductiv¬ 
ity,^ and charge-density-wave and s-wave superconduc¬ 
tivity.!^ Our results expand these examples to the case 
of magnetically ordered systems. 


IV. EXACTLY SOLVABLE LIMIT 

In this section we consider quenches for which the 
nonadiabatic dynamics of the pairing amplitude can 
be found exactly. Specifically, we consider quenches into 



FIG. 2: (color online) Same as Fig. [^for Tc = 0.95Tm 



FIG. 3: (color online) Same as Fig. [^for Tc = 0A5Tm- 


a state with a zero SDW order parameter, rriz = 0. For¬ 
mally, this limit can be realized for the quenches of the 
SDW coupling constant gm —t 0. 

The equations of motion are governed by the following 
Hamiltonian, 

H = 2 ^ ^ (A,AT- + A,A:+) , (29) 

p2 pi 


which we write in terms of the Anderson spins, where 
i = c,f. 


A. Integrability criterion 


In order to demonstrate the exact integrability of 
the model (29), we adopt the method developed by 
Yuzbashyan et al. [65166) . The central role in finding 
the dynamics of the pairing amplitude is played by the 
Lax vector. In order to identify the expression for the 
Lax vector for our problem, we first introduce 


a^ = 



i = c 

i = f 


(30) 
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From Eq. 
(p^ as 


it follows Ai = ajA. Using (30) we rewrite 


2^ — 9sc 'y ^ ®i^pi — 9sc y ' ^pi) (31) 

pi pi 


where we used = —at, and redefined the anomalous 
components of the pseudospins according to ajK^f —>• 
Kpf. Clearly, this transformation leaves the Poisson 
brackets invariant: 




ibc 


k: 


pi- 


(32) 


Thus, the Hamiltonian (29) can now be rewritten as fol¬ 
lows 


H = (33) 


Consider the following Lax vector 

^ sp? ysc 
PJ 


(34) 


The Poisson brackets for the components of C are ob¬ 
tained using (32): 


{C'^{u),C\v)} = 


- C%V) 


(35) 


Due to the commutation relations (35) it follows that 

{C{u),c\v)} = 0. (36) 


This property means that any model Hamiltonian which 
Poisson com mutes with will define an exactly inte- 
grable model!^^^^ Indeed, for the square of the Lax vec¬ 
tor we readily find 



FIG. 4: (color online) Imaginary part of the root of C^{u) = 0 
as a function of spin-density-wave magnetization The 

parameters are Ai = O-OISef and 5k = 0 (no Fermi surface 
mismatch). 


One then finds 


^ _ o \" ^pi I \" ^pi I 4 
^ u - ^p* ^ (m - Cpi)^ 9‘^sc ' 


pz 


pz 

we can write 


From the definition of ' 

9sc ^ ^ ‘^Cpi'^pi “1“ const. 


(40) 


(41) 


Since £^(u) is conserved by evolution due to (35) and 
all Zipi Poisson commute with (u ), it implies that Hpi 
are integrable. Furthermore, since Hm^=o is given by a 


linear combination of Hpi (41), it also commutes with 


C^{u) and therefore is integrable. 


B. Lax roots 


c^(u) = W + A 

(w-Cpi)(M-'Cqj) 5sc 
4 -^pi 

9sc q 7 ~ ^pi 


(37) 


The first term and the last terms should be rewritten as 
follows, 


Kpi ■ Kqj 4 Kp^ 

pi qj (“ ~''P*)(^ ~ ^qj) 9bc U — ^pi 

^2 

^ o'Sp '^P» , ^Pi 

^U-^pi ^ {u-^piY’ 
pz ^ pz ^ / 


(38) 


where 


■ (39) 


qv^p i/i 


Cpi ^qf 9sc 


To determine the value of the pairing amplitude at long 
times, we need to compute the imaginary part of the Lax 
roots governed by the equation 


C^{u) = 0. 


Using (34) we rewrite (42) as follows. 


£^(u) = ±i£x(u), 


(42) 


(43) 


where we took into account that initially all y compo¬ 
nents of the pseudospins are zero, = 0. Furthermore, 
since there is no mismatch between the Fermi surfaces, 


we have ^pj = %— M = —■^p.c = Cp- From Eq. (21) it 


follows Ek = \/^k + where Aq is a supercon¬ 

ducting order parameter to be specified below. Keeping 
in mind that the z-components of Kp\ have not been 
rescaled, for (20) with ^/ = —= C we have 


Ck 




^c(k,t = 0) = ^, Kf(k,t = 0) = -^. (44) 


2E, 


2Ek 
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Similarly, for the anomalous components (221, we find 

An 


K^iKt = 0)=a, + 

iy;(k,f = 0) = «. 


(45) 


Thus, we observe that 

_ ^pi 


= 

P* 2E^ 


K^- = — 
P* 2E„ 


1 - 


m; 

E"^ 

P 


(46) 


a. Self-consistency conditions. For quenches into 
the purely superconducting state, = 0, a different 
equilibrium value of the pairing amplitude Ag is deter¬ 
mined by the BCS self-consistency condition. For a given 
value of the superconducting coupling we have 


gsc p 


1 




(47) 


As a next step, we introduce the function g'scir^z) which 
formally enters as a “new” coupling constant. The equa¬ 
tion which determines this function reads 


5sc(Wz) 


= E 


1 




+ TO? + A? 


(48) 


By comparing (47l with (48) we see that g'^A^) = gsc, so 
that the imaginary part of the Lax root Im[u] = Ag, as 
it should be for equilibrium. 

b. Equation for the Lax roots. Using expressions 
(46) we now rewrite (43) as follows. First, momentum 
summations are replaced with integrals over e = p^/2 —/r 
according to the formula: 

OO 

= f F{e)de 


(49) 


-M 


and vp is the density of states at the Fermi level. For the 
^-component of the Lax vector (34) using (49) we have 


OO 

Fz{u) = -vp / X! 

,, A=± 


Xede 


2{u- Xe)E{e) 


(50) 


= —vpjd — vpu / 


de 


\=± 


2(u — Xe)E{e) ’ 


where we employed Eq. (48) and introduced the param¬ 
eter P, which describes for the magnitude of the quench 


P — ^Pp I 

\gsc g'sc 


1 


1 


(51) 


Note since g^ > 5sc parameter P always remains positive, 
/3 > 0. We can now use expression (50) to rewrite Eq. 



FIG. 5: (color online) Quenched dynamics of the s* pairing 
amplitude as a function of time (5 is a level spacing) for initial 
conditions with a nonzero spin-density-wave magnetization 
roz = SAi (top panel) and = 0.5Ai (bottom panel). 


(43) as 


-/? 

U ± jAg 



iAgm? \ 7 
u±iAo) J ^ 


de 

2{u — Xe)E(e) 


■ (52) 


As we can see from analyzing this equation for ruz = 0, 
there will be only one complex root u = ±*Ag. For 
nonzero ruz we therefore need to find all the complex 
roots of this equation. The imaginary parts of these roots 
will determine the value of the superfluid order parameter 
at long times, see Figs. and 


V. DISCUSSION AND PERSPECTIVES 

In this work we have described nonequilibrium kinet¬ 
ics of order parameters in the context of multicomponent 
superconductors with the emphasis on the iron-pnictide 
systems. We have found, that out of equilibrium, the cou¬ 
pling between competing superconducting and magnetic 
orders occurs not only by virtue of self-consistency con¬ 
ditions but also dynamically. This becomes essentially 
transparent in the pseudospin representation of equations 
of motion for the Green’s function. In particular, as can 


























be seen directly from Eq. (26), precession of the super¬ 
conducting pairing amplitude is strongly affected by the 
dynamics of magnetic order, which then itself back acts 
on the m(t). 

Insofar as our analysis is limited to the collisionless 
regime at time scales satisfying 


TA < t < Tin, (53) 

where ta = 1/A, Ti„ is the time scale of inelastic scat¬ 
tering processes in the collision-dominated regime. The 
latter can be found from the golden rule by passing to the 
Bogoliubov quasiparticle representation where the ma¬ 
trix elements of the transition probabilities in scattering 
are dressed by the coherence factors. Following the early 
works of Eliashber^^ and Galaiko,!^ one estimates 


r-n {gscyFfTljEF, (54) 


which is essentially a Fermi liquid expression for the time 
scale of electron-electron collisions. It is expected that 
the power-law decay of the order parameter crosses over 
to exponential behavior oc exp(—t/Tj„) once the system 
enters into the collision-dominated regime of relaxation. 

In the modeling we have adopted the band model of 
FeSC, which is certainly suitable for clean 122-systems 
such as isovalently P-doped BaFe 2 (Asi_a;Pa;) 2 - It is 
worth pointing out that in a more general formulation 
(for example, within the three-band model) additional 
features may appear, in particular, possibly different 
branches of collective excitations. It is of clear experi¬ 
mental relevance to r evisit the same problem for the dis¬ 
order model of FeSC,l^2l2n! is more appropriate for 


the Co-doped case Ba(Fei_a;Coa;) 2 As 2 . However, phys¬ 
ically perhaps the most interesting question is to study 
the nonequilibrium dynamics near the quantum critical 
point, namely, near the end point of SDW order under 
a superconducting dome. Such quantum criticality was 
revealed fro m the measurements of the London pene¬ 
tration deptt|AlI21 and attributed to the fluctuations of 
SDW order at the onset of the transition into the coexis¬ 
tence phase.l^ The dynamics of magnetization near such 
a quantum critical point has been recently addressed in 
the framework of time-dependent Ginzburg-Landau the¬ 
ory,!^ however, the nonadiabatic regime still needs to 
be systematically investigated. In general, post quench 
prethermalization at a quantum critical point may ex¬ 
hibit nontrivial dynamical scaling .1^ The case of iron 
pnictides is very specific since the magnetic quantum crit¬ 
ical point is surrounded by a superconducting state with 
gapped quasiparticles, and consequently scaling of the 
response functions may be governed by entirely differ¬ 
ent dynamical exponents. Finally, one should seriously 
look at the role of degrees of freedom associated with the 
Ising nematic order parameter!^ that was left behind in 
our picture. All these questions will pave the way for 
future research in this field. 
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